Focusing of laser-generated ion beams by a plasma cylinder: similarity theory and the 

thick lens formula 
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It is shown that plasma-based optics can be used to guide and focus highly divergent laser- 
generated ion beams. A hollow cylinder is considered, which initially contains a hot electron pop- 
ulation. Plasma streaming toward the cylinder axis maintains a focusing electrostatic field due to 
the positive radial pressure gradient. The cylinder works as thick lens, whose parameters are ob- 
tained from similarity theory for freely expanding plasma in cylindrical geometry. Because the lens 
parameters are energy dependent, the lens focuses a selected energy range of ions and works as a 
monochromator. Because the focusing is due to the quasineutral part of the expanding plasma, the 
lens parameters depend on the hot electron temperature Te only, and not their density. 

PACS numbers: 



I. INTRODUCTION 

Laser-driven ion sources P, S S S IE IE consid- 
ered to be the hot candidates for various important ap- 
plications in nuclear physics, medicine, biology, material 
sciences, plasma field tomography |EIE^3- When multi- 
terawatt laser pulses are shot on solid state targets, copi- 
ous amounts of multi-MeV ions - both protons and highly 
charged heavier ions - are generated These laser- 

generated ion beams have picosecond durations and orig- 
inate from a few micrometer wide virtual source. How- 
ever, the laser-generated ions are highly divergent and 
usually are emitted within a cone with some 10-30 de- 
grees opening angle. In addition, they have broad energy 
spectra. These facts may impede numerous applications 
for the laser-generated ion beams unless appropriate op- 
tics and monochromatizing systems are developed. 

Because of their high divergence, one needs very strong 
fields to collimate the ion beams. Such fields exist only 
in plasma. However, one cannot exploit the standard 
technique of self-induced magnetic plasma lensing that 
is widely used to focus conventionally accelerated ion 
beams. The reason is that the laser-produced ion beams 
are charge neutral, i.e. they contain electrons that com- 
pensate the ion charge and current. 



II. SIMILARITY THEORY OF EXPANDING 
PLASMA 

In the present work we consider ion beam focusing by 
plasma which already contains a quasistatic electric field. 
The experimental configuration is the following. Fig. ^ 
[l^: a laser produced ion beam originates at the point 
z = ~D, X = y = and propagates toward a hollow 
metal cylinder of the radius R and length L, L > R. 
The axis of the cylinder coincides with the z-axis. At the 
same time, the second laser pulse is shot at the cylin- 
der. This second pulse generates a population of hot 
electrons, which penetrate through the metal and spread 



FIG. 1: (color online). Geometry of the ion optics element 
for focusing of laser-generated ion beams. One laser beam 
generates an ion beam from the rear side of the irradiated 
foil. Another laser beam hits a hollow cylinder, where a hot 
electron population is generated. The cylinder has the radius 
R, the length L and is located at the distance D from the 
foil. The radial electric field of the plasma inside the cylinder 
coUimates the ion beam. 



very fast over the inner surface of the cylinder. They exit 
into vacuum and generate a cloud of space charge at the 
inner surface. The electric field of this space charge is 
large enough to ionize the material and to create plasma. 
As a result at the initial moment we have a cylindri- 
cal plasma layer with high electron temperature Tg and 
low ion temperature « 0. The plasma begins to ex- 
pand toward the cylinder axis due to the TNSA (tar- 
get normal sheath acceleration) mechanism Nor- 
mally, the cylinder surface is covered by a thin layer of 
hydrogen-rich substances. Being the lightest ions, pro- 
tons are accelerated first and the plasma is usually an 
electron-proton one regardless of the particular chemical 
mixture of the cylinder itself. 

Plasma dynamics is described by the couple of Vlasov's 
equations for electrons and ions and the Maxwell equa- 
tions 
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where fe and /i are the electron and ion distribution 
functions respectively. We do not include magnetic field 
in 1^ since the cylindrical symmetry of the expansion 
prohibits magnetic fields generation. 
The initial conditions for Eqs. ^ are 



/e(t = 0,p,r) = fo{cp/T,,r/R,r/d), 
Mt = 0,p,r) = Fo(r/i?,r/d)5(p), 



(2) 



where /o and Fq are initial distributions of electrons and 
ions, d and R are cylinder thickness and radius respec- 
tively, R > d. For the initial distributions one reads 

//„(cp/r.,r/fl.,/<i)<ip<iR = 2,«.L„.. 

J Fo{r/R,r/d)S{p)dpdR = 2TTRdLne/Z^, (3) 

where the cylinder length L ^ R, is the average elec- 
tron density. The multiplier 1/Zi is due to the condition 
of plasma charge neutrality. 

In the following we consider the case of relativistic 
electron temperatures and non-relativistic ions: MiC^ 3> 
Te ^ iTLeC^ ■ Thus, wc cxploit the ultrarelativistic ap- 
proximation V = cp/|p| for the electron velocities, while 
for the ions we get v — p/Mi. Here Mi and rrie are the 
ion and electron masses respectively. 

We introduce new dimensionless variables 



f = r/i?, i^t/T, T^R/cs, E = ei?E/Te, (4) 



where Cs — y^Te/Mt is the ion sound velocity. Since 
the electrons are ultra-relativistic and ions are non- 
relativistic, we are forced to introduce different dimen- 
sionless normalization for electron and ion components: 

Pe = CPe/Te, = [Tjcf fjn„ (5) 

Pi = pi/^Iufe, /, = {M,T,f'^ fjn,. (6) 
We rewrite the Vlasov equations in these variables: 



dfe{i,Pe,r) Pe 5/e(t,Pe,f) 9/^ (i, Pe, f ) „ 
Ctr h T-TT E — = 



dt 



dv 



dpe 



d Mi, Pi, r) ^ df,{i. Pi, r) ^ ^f^ {i, pi , f ) 

■ + Pi — hAE = U, 



di ' ^' Or dpi 
az) V • E = 47re ( Z, [ f, d^p - f fe d^p 



(7) 



The normalized Vlasov-Maxwell equations lO reveal that 
the the system dynamics depend on five dimensionless 
parameters. The first parameter is the ion charge Zi. 
The next two parameters are the normalized sound speed 
Oc = Cs/c and the normalized Debye length ao — 
A|,/47ri?^, where Xj^ — AirTe/ e^Ue- These two param- 
eters define plasma dynamic properties. The remaining 
two parameters d/ R and L/ R come from the initial sys- 
tem geometry. We are interested in the cylindrical geom- 
etry and drop out the parameter LjR^ -|-oo. 

Thus, the parametric dependencies can be written as: 



/e = 



t r cp d 

t'R'T^'R' 
t r 



Zi,ac, ao 
p d 



(8) 
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where and fi are universal functions. Eqs. |(SJ) already 
can be used to state exact scaling laws. The requirements 
Q!c — const, an = const and d/R — const do not fix 
all the dimensional parameters of the problem, and this 
allows to scale experimental results. 

Yet, the most interesting scalings are obtained in the 
limit Qfc ^ 1 and ajj ^ I. Assuming ac — > one obtains 



Pe 
|Pe 



dfe{i,Pe,r) 
di- 



E.^MlP^^O. (9) 

OPe 



This means that the electron distribution fimction can 
be written as 



fe = Fe (|Pe| - 4>, t, Zi,d/R, Ud 



(10) 



where Fe is a universal function. Eq. (|10|l means that 
the electron fluid has the same effective temperature at 
all points. 

The formal limit — > coincides with the quasineu- 
trality condition 



(11) 



Z^ / f^d•'p= / fed'p) 



Since an is a factor in front of the highest derivative in 
(UI), the quasinetrality condition Hill) is violated within 
the narrow Debye sheath layer of the width cx Ad. Be- 
ing very important for problems like ion acceleration this 
area hardly plays any role in the ion focusing. Because 
of its narrowness, only the small amount of ion beam 
on the order of Xd/R ^ 1 would be influenced 

by its flelds at any particular moment. We neglect this 
influence. 



III. ION FOCUSING BY HOLLOW PLASMA 
CYLINDER 

In order to describe the focusing, we study properties 
of the quasineutral part of the expanding plasma cloud. 



3 



The quasineutrality Ue ~ ZiUi is guaranteed as long as 
ajj —^ 0. Here rtg and are the electron and ion den- 
sities, Zi is the ion charge state. At the same time, the 
plasma density and consequently the electron pressure 
Pe — UeTe vary along the cylinder radius. The pressure 
gradient is counterbalanced by the radial electric field 



E 



1 



-VP. 



(12) 



which is developed inside the plasma to satisfy the 
quasineutrality Hll() . Because the electron pressure gra- 
dient is directed off axis, the developed electric field is 
directed toward the cylinder axis. It is this field that 
focuses the injected ions. 

Because of the cylindrical symmetry, we neglect any 
dependencies on the azimuthal angle on the longitudinal 
coordinate z within the plasma. Then, all distributions 
depend on the radius p = \J + 'tp- only. To obtain 
a closed system of equations we take into account the 
energy conservation law 

^TeA^e + TTML ^ v^/. (i, V, p) dv pdp = ^T^ (O)iVo, (13) 

where iVg is the number of hot electrons and 
(3/2)re(0)A^e is the laser energy absorbed in the cylin- 
der and stored in the hot electrons. Eq. ifT^ neglects 
the energy accumulated in the electromagnetic plasma 
fields. This assumption is correct provided that the De- 
bye length is much smaller than the cylinder radius P, 
i.e. for ^ 1. Eqs. @ and H1U|I show that the electron 
temperature is equal at all points of the plasma. Thus, 
the energy conservation law (|13|l is sufficient to describe 
the electron dynamics. 

The initial ion distribution is 



h{t = 0, p, p) = 27ra,<5(p)Po(p/i?, diR)- (14) 

where Gi is the initial surface density of ions participating 
in the plasma expansion. Because of the quasineutrality 
condition pi(l . we have iVe = l-KRLZiUi. 

We introduce the dimensionless time-dependent elec- 
tron temperature Tit) = Te{t)/Te{0) and the ion velocity 

V = V/Cs- 

The ion Vlasov equation and Eqs. H12|) - (|13|l rewritten 
in the dimensionless variables take the form: 



Eqs. (|15() - (|18|l contain no dimensional parameters what- 
soever. As a consequence, the functions T, / and E are 
universal, i.e., they are not affected by specific values 
of d, P, L, Ui and Te(0). This gives us an opportunity 
to develop a meaningful similarity theory describing the 
guidance of laser produced ion beams. 

From the normalizations we conclude that the elec- 
tric field E developed in the plasma is 



E 



TejO) 

eP 



E{t/T, p/R, d/R). 



(19) 



where E is a universal function. It does not depend on 
the plasma density, but is determined by the hot electron 
temperature and the cylinder geometry only. This result 
is valid as long as the Debye length is much smaller than 
P. This means that the uncompensated charge density 



e5n = eiZiUi 



^-^m/r,p/R,d/R) 



VE 

47r 



47reP2 



(20) 

is much smaller than the electron density. 

When the laser produced ion beam enters plasma in- 
side the cylinder, it is defiected by the electric field H19|l . 
We suppose that the beam has a lower density than the 
plasma inside the cylinder and thus the beam own fields 
can be neglected. 

To describe the beam ion guiding in plasma we consider 
ions with the charge state Zh, mass and the initial 
energy £b being focused by the potential 



(fi = —nr'^eSTiQ, Sn^ — STi{t/T, p — 0, d/ R) (21) 

Notice that the charge density Suq depends on time. 
However, for the most interesting and important case 
the beam ions pass the cylinder plasma during the time 
L/ub <C T, where Ub ~ \J 2£b /Mb and r — R/cs is the 
plasma evolution time. In this case the dependence of 5n 
from time t can be neglected. 

Now we are able to estimate infiuence of the non- 
neutral Debye sheath with the width Ad on the beam 
ions motion. This area propagates with the velocity cx Cs 
and carries the electric field P„g cx y/nXTl. The radial 
momentum of a beam ion is changed by the value 



di 



E 



-E 



Vln rie 
0, 



w^fS,^,p)d^df - 3(1 -T), 



with the initial condition 



/,(t = 0,v,p) = Po(p,d/P)(5(v). 



(15) 
(16) 

(17) 
(18) 



Ap""^ cx eZiy^UeTe cx AEnqMb, (22) 

where ASnq oc j^Te. 

The change of a beam ion radial momentum due to the 
interaction with the quasineutral part of plasma is esti- 
mated as 



ApJJ_ cx 



R Ub 



cx 



V^^gMb, 



(23) 
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where cx (^) ' f T,. 

From Eqs. ^ and ^ one sees that Ap"^ » Ap]_. 
Therefore the ions passing through the non-quasineutral 
edge are stronger deviated than those interacting only 
with the quasineutral plasma region. However, because 
the Debye sheath is narrow, the relative number of the 
strongly declined ions is small and these ions are deviated 
to different points of space. For these reasons the Debye 
sheath at the edge of the expanding plasma does not 
contribute to the ion focusing. It scatters the beam ion 
instead. 



Note that when the electrons pull the ions behind them- 
selves (positive pressure gradient) there is an electron ex- 
cess at the axis z, i.e., Suq < 0. The constant C in H29|l 
is obtained from the continuity conditions of the action 
S at the front and rear sides of the plasma cylinder. 
Thus, we arrive at the thick lens formula 

{Z-g){Z' + h) = -f (30) 

where g = — (L cose)/ (e sine), h = g + L, f = L/{esme) 
and 



IV. ION LENS FORMULA 

To investigate focusing properties of the potential (p we 
use the well known analogy between the geometrical op- 
tics and the classical mechanics. The optical length corre- 
sponds to the action S in the Hamilton-Jacoby equation 



spo r 

M 



dtS + H{VS,r) = 0. 



(24) 



If ions in vacuum are injected at the point x ^ y = 0, 
z ~ Z then the S function in vacuum is 



y/ZbSh{t/T, 0, d/R) L 

2 UfcT 



(31) 



In our derivation of Eq. (|30|l we neglected the change of 
plasma parameters during the time the beam ions need 
to pass the cylinder. This means that our analysis is valid 
if e< 1. 

A parallel beam of ions is obtained if Z' = oo. This 
condition is satisfied for e <C 1 if 



D 



L 
72- 



(32) 



S = -£bt + y/2Mb£b{{z-Z)^ + p^) 
« -£bt + ^2Mb£b{z -Z) + 

z ~ Z 



(25) 



Mb£b 



In our geometry, Z = —D. 

If ions are focused at the point a; = y = 0, z = Z' va 
vacuum, then the action S is 



S = const - £bt - ^y2Mb£b{{z - Z')^ + p^) 



(26) 



ist - £bt + ^/2Mb£b{z - Z') 



Mb£b 



The beam ion motion inside the plasma cylinder is de- 
scribed by the Hamiltonian 



H 



2Mb 



eZbifi. 



(27) 



The solution of the Hamilton-Jacoby equation inside the 
plasma can be expanded as 



Thus, the plasma element coUimates ions with the energy 



£h oc ZhTf 



LD 



(33) 



It is worth mentioning that the energy of the coUimated 
ions strongly depends on the dimcnsionless parameter 
LD/R^ and can be significantly larger than the initial 
electron temperature. 

It is easy to see that for the ion focusing be practical 
the electron temperature has to be of the order of 
several Me Vs. Such electron temperatures are routinely 
produced by multi-terawatt lasers. 

Would be the plasma inside the cylinder stationary, 
then only ions with the selected energy (|33|l are coUi- 
mated. However, the plasma is non-stationary with the 
characteristic evolution time r = R/cg- The relative 
change of the plasma parameters during the ion passage 
time through the cylinder is of the order of L/vbT <^ 1. 
This small parameter defines finally the finite energy 
spectrum width S£b of the focused ions: 



1 



S = const - £bt + y/2£bMbZ + -I3{z)p^ + (28) 
where the function (3 is 



f3{z) = \/ —2TrZbMbe^Sno tan 



(29) 



S£b L 
— cx cx 

£b VbT 



(34) 



It follows from H34|l that the plasma cylinder works as a 
good monochromator if D 3> i. 

To avoid any confusion we emphasize that Eq. I|34(l 
describes the quality of a small aperture ion beam only. 
Of course, different parts of the lens coUimate ions of 
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different energies. For large aperture ion beams the en- 
ergy spectrum width of the focused ions wiU be large, 
A5h/£f, cx 1. 

In the preceding part of the paper we consider the fo- 
cusing by the area p ~ 0. It is important that the theory 
can easily be generalised for the focusing by the cylindric 
area near p = po < R. To do so we introduce a new 
potential 



Since the plasma lens density gradient can be very large 
this condition can be much weaker than ni, <ti Ui . 



CONCLUSIONS 



•idrir 



TejO) 



E{t/T,po/R, d/R). 

(35) 

According to Eq. (fH^ the potential (fro gives the right 
value of the electric field at p = pq. Thus the focusing 
by p = po is obtained from Eq. H30() with e ^ 1 by the 
substitution Suq —f Suro ■ 

Until now we have assumed that the density of the ion 
beam focused is so small that it does not affect the fo- 
cusing field of the lens. To find the validity condition for 
this approach we have to consider the propagating of an 
ion beam with a given density profile nfc(p) through the 
lens. Using the quasuneutrality condition for the system 
" the lens plasma -I- the beam" one can easily find the the 
focusing by the area around p — po is not disturbed by 
the ion beam if 



In conclusion, we have developed a closed similarity 
theory of a hollow cylinder as a plasma element for ion 
beam guiding. Significantly, the beam ions are focused by 
the quasinuetral part of expanding plasma rather than by 
strong electric fields in the non-quasineutral leading edge 
of the expanding plasma cloud. The thick lens formula 
has been obtained with explicit scalings for all of the 
parameters. We show that the plasma lens collimates 
only ions with a quite definite energy and may be used 
for monochromatization of the laser-produced ion beams. 
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